Reflectance Spectroscopy as a Novel Tool for Thickness Measurements of Paint Layers

A major challenge in heritage science is the non-invasive cross-sectional analysis of paintings. When low-energy probes are used, the presence of opaque media can significantly hinder the penetration of incident radiation, as well as the collection of the backscattered signal. Currently, no technique is capable of uniquely and noninvasively measuring the micrometric thickness of heterogeneous materials, such as pictorial layers, for any painting material. The aim of this work was to explore the possibility of extracting stratigraphic information from reflectance spectra obtained by diffuse reflectance spectroscopy (DRS). We tested the proposed approach on single layers of ten pure acrylic paints. The chemical composition of each paint was first characterised by micro-Raman and laser-induced breakdown spectroscopies. The spectral behaviour was analysed by both Fibre Optics Reflectance Spectroscopy (FORS) and Vis-NIR multispectral reflectance imaging. We showed that there is a clear correlation between the spectral response of acrylic paint layers and their micrometric thickness, which was previously measured by Optical Coherence Tomography (OCT). Based on significant spectral features, exponential functions of reflectance vs. thickness were obtained for each paint, which can be used as calibration curves for thickness measurements. To the best of our knowledge, similar approaches for cross-sectional measurements of paint layers have never been tested.


Introduction
In recent decades, a wide variety of scientific techniques have been tested and optimized for the study of cultural heritage (CH) objects. The need to preserve the material integrity of works of art has directed the research toward defining non-invasive analytical approaches based on the combined application of methodologies that do not involve sampling or risk of damage to the object. A major challenge in heritage science is the non-invasive cross-sectional analysis of the pictorial stratigraphy in paintings. Thickness measurement of the micrometric layers is essential, for instance, to monitor the removal of surface materials during the cleaning operation [1] or to assess the compactness and adhesion between layers. In a non-invasive approach, when low-energy probes are used, the presence of opaque media can significantly hinder the penetration of the incident radiation, as well as the signal collection from within the examined materials.
One of the most widely used techniques for non-invasive stratigraphic measurements on paintings is Optical Coherence Tomography (OCT) [2,3]. Primarily applied in the field of ophthalmology, OCT is an interferometric method based on a Michelson interferometer, yielding 2-or 3-dimensional tomographic imaging that allows the visualization of the internal structure of pictorial layers with an axial resolution ranging from 1 to 10 µm (in air). The incident radiation is backscattered by the material and the optical interference is observed whenever the signal superposes with the reference beam, within the coherence length of the light source. The measurement is based on the detection of signals generated at the interfaces between different media-i.e., when the incident radiation experiences a refractive index (n) mismatch. OCT has proven particularly effective in probing materials that are semi-transparent in the near-infrared (NIR) spectral range. This includes most varnishes applied by artists on the painting surface with a protective and/or aesthetic function [4]. By combining an OCT setup with confocal microscope optics, which enables the beam focussing inside the material rather than on the outer surface, even highly reflecting coatings can be measured [5]. Pictorial layers, however, are often composed of pigments with dispersion and/or absorption properties that do not allow their thickness to be assessed by OCT. Moreover, the n-mismatch causes a delay in the optical path of the reference beam and, therefore, the optically measured distances must be corrected to geometrical distances by dividing them by the refractive index of the material. In the NIR, the refractive index of semi-transparent materials used in paintings is conventionally given as 1.5. However, pictorial layers (pigment dispersed in the binder) are often highly heterogeneous and exhibit variable optical properties that result in different n values. Therefore, if n is not known, the correct thickness of the painting layers is not achievable.
In recent decades, alternative methods to OCT have been proposed for the noninvasive in-depth analysis of paintings. Among others, Terahertz imaging [6,7] has proven effective in yielding 3D data sets of interfaces and projections of paintings in the presence of absorbing species. However, the low axial resolution achievable with THz radiation makes this method unfit for micrometric measurements of pictorial layers.
The use of Nuclear Magnetic Resonance (NMR) to obtain stratigraphic information on easel and wall paintings is also well-documented in the literature [8,9]. NMR profiling was tested to investigate both signal intensity and transverse relaxation time distribution, showing that the dependence of signal intensity on relaxation times makes the interpretation of the stratigraphic information difficult [10]. The NMR-sensitive volume averages the effect of irregularities in the layers and the signal from adjacent layers. In addition, the application of this method is hampered by the lack of application-specific operating software, while the low mass sensitivity resulting from low NMR frequencies results in long measurement times.
More recently, Nonlinear Optical (NLO) techniques [11] have been successfully used for cross-sectional analysis of a wide variety of artistic materials, including paint and varnish samples. The combined application of different NLO modalities allows for the acquisition of compositional and structural information based on the detection of fluorophores (by Multi-Photon Excitation Fluorescence, MPEF), crystalline or highly organized structures without inversion symmetry (by Second Harmonic Generation, SHG), or local differences in refractive index, i.e., interfaces (by Third Harmonic Generation, THG). While for MPEF and SHG, the main limitation in the stratigraphic analysis of paintings is the presence of highly diffusing and/or absorbing media (pigments) [12], the applicability of THG is confined to layers of transparent material, forward detection being the only possible configuration [13].
A cutting-edge methodology recently proposed for cross-sectional analysis in paintings is photoacoustics, which in a sense, can be considered complementary to OCT, as it takes advantage of the presence of non-transparent materials [14,15]. Acoustic waves are generated by the absorption of the radiation emitted by an intensity-modulated pulsed laser. The exponential attenuation of acoustic waves in the frequency domain, which depends on the absorption coefficient of the medium and the propagation path, can be exploited to measure the thickness of the examined material. Although early applications reveal the technique's potential [16], to date, its use is limited to specific cases only.
Given the above, it can be stated that, at present, there is no technique that can uniquely and non-invasively measure the thickness of heterogeneous and optically opaque materials such as pictorial layers.
In this work, we explored the feasibility of achieving stratigraphic information of painting layers from their reflectance spectra measured by Diffuse Reflectance Spectroscopy (DRS) [17]. In heritage science, DRS is typically applied for the analysis of paintings in multiand hyper-spectral imaging modes or using fibre optics for point-wise measurements. The main objective is typically the identification and mapping of pigments and binders based on the absorption properties of electronic and vibrational transitions of molecules [18]. In the imaging mode, the use of the NIR spectral range allows for the visualization of hidden details underneath the painted surface related to the artistic working process, such as underdrawings and underpaintings. Diffuse reflectance is defined as the ratio of the irradiance of light reflected back to the detector to the irradiance on the surface of the object, as a function of wavelength. The measured light backscattered from the object includes contributions from the air/surface interface and varies with illumination and collection geometry.
The spectral reflectance behaviour of pigments and pictorial layers has been extensively studied [19][20][21]. It has been shown that the reflectance signal measured from pigment mixtures in paintings is the result of the nonlinear combination of the reflectance of the individual pigments. To cope with the complexity of spectral data interpretation, as well as to reduce the high dimensionality of DRS imaging datasets, new approaches based on artificial intelligence (AI), e.g., deep neural networks (DNN), have been recently explored [22].
In this preliminary study, we tested the proposed approach by examining single layers of pure paint, thus avoiding the use of optical models to predict the reflectance of pigments in mixtures [23] and not taking into account the influence of the surface roughness on the spectra [24]. A mock-up was created for this specific purpose: ten acrylic paints were laid with increasing thicknesses, ranging from 50 to 350 µm, on both white and black backgrounds. The chemical composition of each paint was first characterised by micro-Raman [25] and laser-induced breakdown spectroscopy (LIBS) [26]. The thickness of each layer was then measured by OCT, taking as a reference the portion of the visible substrate at the edge of the paint layer. The spectral behaviour was analysed by both Fibre Optics Reflectance Spectroscopy (FORS) and Vis-NIR multi-spectral reflectance imaging. Finally, for each acrylic paint, meaningful spectral features were identified to assess the dependence of the reflectance on the layer thickness, thus obtaining non-linear fitting curves.
To the best of our knowledge, no similar approaches have been explored before for cross-sectional measurements of paint layers.

Chemical Characterization of the Acrylic Paints with LIBS and Micro-Raman Spectroscopy
The chemical composition of the ten acrylic paints was assessed by LIBS and micro-Raman spectroscopies (Table 1), using the information reported in the literature [27][28][29][30][31][32] and NIST [33] and IRUG [34] databases. Further information on the chemical composition of the phthalocyanine paints (PBC, PBL and PGL) can be found in our previous work [35]. Figure 1 shows the LIBS (top) and micro-Raman (bottom) spectra of cobalt blue (CB), cadmium red (CR), cadmium yellow (CY) and primary blue cyan (PBC) acrylic paints. The spectra measured on the other analysed paints are displayed in Figure S1 in the Supplementary Material. For all paints, the chemical composition (Table 1) agrees with what was declared by the manufacturer (Table 2), except for the absence of titanium dioxide in permanent green light (PGL) and the presence of additional components in most of the analysed paints, ascribed to the binder and fillers. Specifically, atomic emissions of Mg, Si, Ca, Al, Sr, Ba and Na detected by LIBS are ascribed to fillers such as kaolin (Al 2 Si 2 O 5 (OH) 4 ), gypsum (CaSO 4 · 1 /2H 2 O), carbonates (CaCO 3 , MgCO 3 ), glass powder and barite (BaSO 4 ) [31]. The molecular bands of CN (Violet band), CH and C 2 (Swan bands) are due to the acrylic binder, as well as to the organic pigments. The LIBS results obtained in these findings agree with the ones by micro-Raman spectroscopy. Bands from calcium sulfate (1007 cm −1 ), calcium carbonate (1085 cm −1 ) and barium sulfate (985 cm −1 ) are observed in each acrylic paint and attributed to the fillers. Additional bands at 482, 600, 620, 837, 841, 1106, 1150-1200, 1240, 1305, 1449, 1452, 1728, 2411 and 2800-3100 cm −1 are attributed to the constituents of the polymeric binder [27][28][29][30][31][32]. Table 1. Summary of elemental and molecular composition of the acrylic paints as found by LIBS and micro-Raman spectroscopies. The main elemental components and the Raman characteristic bands are indicated in bold. LIBS results reported in Table 1 show the presence of the main paintings markers such as Co for cobalt blue (CB), Cd for cadmium red (CR) and cadmium yellow (CY), Cu and Ti for permanent blue light (PBL), Cu for permanent green light (PBC and PGL), CN and C 2 for primary red magenta (PRM), Ti, CN and C 2 for primary yellow (PI), Ti for titanium white (TW), and Zn for zinc white (ZW).

Thickness Measurements with OCT
Four xz tomograms (8 × 0.6 mm 2 , pixel size 3.5 µm 2 ) were acquired in each painted area. Given the low transparency of most of the analysed paints, the layer thickness was measured by taking the signal generated at the interface air-background (visible at the edges of each painted area) as a reference, as shown in Figure 2. The thickness of each pictorial layer was calculated as the average over 20 values, resulting in a range between a minimum of 45 µm (area 1) to a maximum of 350 µm (area 5) for all acrylics. For each thickness, the error, i.e., the standard deviation, resulted below 6 µm for all areas, demonstrating the micrometric homogeneity of the paint layers. Only three acrylics, namely PBC ( Figure 3), PRM, and PY, showed sufficient transparency to enable the evaluation of their refractive index, which was calculated by dividing the thickness measured with OCT by the real one. The resulting n values at 1300 nm, i.e., at the OCT radiation wavelength, are in the range of 1. 35

Thickness Measurements with OCT
Four xz tomograms (8  0.6 mm 2 , pixel size 3.5 µm 2 ) were acquired in each painted area. Given the low transparency of most of the analysed paints, the layer thickness was measured by taking the signal generated at the interface air-background (visible at the edges of each painted area) as a reference, as shown in Figure 2. The thickness of each pictorial layer was calculated as the average over 20 values, resulting in a range between a minimum of 45 µm (area 1) to a maximum of 350 µm (area 5) for all acrylics. For each thickness, the error, i.e., the standard deviation, resulted below 6 µm for all areas, demonstrating the micrometric homogeneity of the paint layers. Only three acrylics, namely PBC (     (1)(2)(3)(4)(5), thickness values are reported in red and calculated as the geometrical distance between the air-paint and the paint-background interfaces, with the la er highlighted by the light-blue line.

Thickness Measurements with Reflectance Spectroscopy
FORS and multi-spectral data were compared for each paint, as shown in Figure 3a-d. In the graphs, the average spectra of each painted area (1→5) are plotted together with the spectrum of the underlying substrate (white or black background). First and second derivatives were computed for all spectra to facilitate the identification of the spectral feature best representing the reflectance dependence on the material thickness. For all paints, the trend of the reflectance as a function of the thickness is expressed by an exponential function, following the equation: where y = R%, y 0 = offset, A = initial value, R 0 = growth constant, and x = layer thickness.
In the case of CB paint, the maximum reflectance R% values in the 808-811 and 710-760 nm ranges were selected for the white and the black background series, respectively, and plotted as a function of the thickness (Figure 4e,f). We noticed that the presence of the black background affects the position of the point of maximum reflectance, causing a blue shift as the thickness of the paint layer decreases and becomes gradually more transparent. In the presence of the white background, however, the point of maximum R remains around 810 nm regardless of the thickness of the paint layer. The resulting exponential fit curves (coefficient of determination R 2 > 0.98) show a good match between FORS and reflectance scanning results.  The spectral feature selected for the analysis in the FORS spectra was not always identifiable in the spectra from the spectral cube due to the significantly lower spectral resolution of the multi-spectral scanner. Therefore, in order to evaluate the applicability of the proposed method in multi-spectral imaging mode, matching key points were found in the two datasets. With this aim, the results of PBC laid on the white background are shown in Figure 5 as an example. The spectral region between 600 and 1200 nm was chosen as significant for our computation: the multi-peak FORS spectra were fitted with a 5th-degree polynomial (Figure 5a) to reconstruct the shape of the multi-spectral spectrum. Maximum reflectance at 950 nm was then considered for both datasets. The resulting exponential fitting functions of the two DRS data show good accordance (Figure 5e). The same maximum was considered for the paint laid on the black background (Figure 5b,d,f) without fitting the FORS spectra to retrieve the same spectral feature. In this case, the reflectance measured on the thickest layer (PBC5 = 284 ± 10 µm) was excluded from the exponential fitting calculation (Figure 5f), since it clearly deviated from the increasing trend, being lower than that of PBC4. This measurable thickness threshold has also been found in other pigments for thicknesses exceeding 270 microns. Remarkably, this limit of detectability was exclusively found in the FORS spectra. This is possibly due to the different measurement configurations between the two DRS modalities, which results in a greater homogeneity of illumination and, therefore, depth of detection achievable with the multi-spectral scanner than with fibre optics.  Optimal agreement between FORS and multispectral data was found in acrylics showing high transparency in the NIR range. As an example, results on PRM on a black background are shown in Figure 6a-c. As for highly scattering pigments, such as ZW shown in Figure 6d-f, an exponential fit curve could be derived only in the presence of the black background. Additionally, in this case, the detection limit found with FORS is Optimal agreement between FORS and multispectral data was found in acrylics showing high transparency in the NIR range. As an example, results on PRM on a black background are shown in Figure 6a-c. As for highly scattering pigments, such as ZW shown in Figure 6d-f, an exponential fit curve could be derived only in the presence of the black background. Additionally, in this case, the detection limit found with FORS is 250 microns for both PRM and ZW, corresponding to the thickness of area 5. As expected, in the presence of the white background, the reflectance of the identified spectral feature was found to decrease exponentially with increasing paint thickness for all acrylics. In contrast, in the presence of the black background, the reflectance was found to increase exponentially. The fitting parameters of function (1) for the ten acrylic paints, together with the considered spectral feature, are summarized in Table 3. In few cases, the high scattering or absorption properties of some pigments made it impossible to identify the spectral feature of interest for the application of the exponential fit. Specifically, the proposed method could not be applied to CR, CY, PBL, TW and ZW laid on a white background and PGL laid on a black background. In all the other cases, the quality of the fitting (coefficient of determination R 2 ) resulted in a value higher than 0.98. Table 3. Fit parameters of the exponential function (1) describing for each paint the dependence of the reflectance of the spectral feature identified with FORS and multi-spectral scanning on the layer thickness.  As expected, in the presence of the white background, the reflectance of the identified spectral feature was found to decrease exponentially with increasing paint thickness for all acrylics. In contrast, in the presence of the black background, the reflectance was found to increase exponentially. The fitting parameters of function (1) for the ten acrylic paints, together with the considered spectral feature, are summarized in Table 3. In few cases, the high scattering or absorption properties of some pigments made it impossible to identify the spectral feature of interest for the application of the exponential fit. Specifically, the proposed method could not be applied to CR, CY, PBL, TW and ZW laid on a white background and PGL laid on a black background. In all the other cases, the quality of the fitting (coefficient of determination R 2 ) resulted in a value higher than 0.98.

Paint
Finally, the proposed method was tested on a sample of cobalt blue laid on canvas (Figure 7a). In the FORS spectrum (Figure 7b), obtained from the average of nine spectra, the maximum at 760 nm, previously identified on the mock-up for CB, resulted in R% = 82 ± 1. The corresponding thickness obtained using the exponential functions derived from both FORS and scanner analysis ranges from 145 to 178 µm (Figure 7c), which includes the actual thickness range measured with OCT, i.e., 150-165 µm (Figure 7d). Finally, the proposed method was tested on a sample of cobalt blue laid on canvas (Figure 7a). In the FORS spectrum (Figure 7b), obtained from the average of nine spectra, the maximum at 760 nm, previously identified on the mock-up for CB, resulted in R% = 82 ± 1. The corresponding thickness obtained using the exponential functions derived from both FORS and scanner analysis ranges from 145 to 178 µm (Figure 7c), which includes the actual thickness range measured with OCT, i.e., 150-165 µm (Figure 7d).

The Painting Mock-Up
The acrylic paints (extra-fine acrylic colours-Maimeri Brera TM , Milan, Italy) selected for this study are reported in Table 2, with their commercial code and chemical composition (as declared by the manufacturer). The surface of the wooden support, a tablet with size 12.5 × 26 × 1 cm 3 , was prepared with a layer of acrylic primer (Lefranc Bourgeois TM , Paris, France), which was accurately smoothed with fine-grit sandpaper. Half of the surface was covered with black acrylic paint (Carbon Black-PBk7-77266), as shown in Figure 8. Each acrylic colour was laid on five adjacent areas of 1 cm 2 (1-5 in Figure 8b) with increasing thicknesses. To this aim, a nearly 50 µm thick tape (Kapton Cypress-DuPont TM , Wilmington, DE, USA) was used as a reference, i.e., depending on the desired thickness, 1 to 5 layers of tape were overlapped on both sides of the area to be painted. A glass spatula was used to spread the pigments in the various thicknesses in order to make their surface as even as possible.

The Painting Mock-Up
The acrylic paints (extra-fine acrylic colours-Maimeri Brera TM , Milan, Italy) selected for this study are reported in Table 2, with their commercial code and chemical composition (as declared by the manufacturer). The surface of the wooden support, a tablet with size 12.5 × 26 × 1 cm 3 , was prepared with a layer of acrylic primer (Lefranc Bourgeois TM , Paris, France), which was accurately smoothed with fine-grit sandpaper. Half of the surface was covered with black acrylic paint (Carbon Black-PBk7-77266), as shown in Figure 8. Each acrylic colour was laid on five adjacent areas of 1 cm 2 (1-5 in Figure 8b) with increasing thicknesses. To this aim, a nearly 50 µm thick tape (Kapton Cypress-DuPont TM , Wilmington, DE, USA) was used as a reference, i.e., depending on the desired thickness, 1 to 5 layers of tape were overlapped on both sides of the area to be painted. A glass spatula was used to spread the pigments in the various thicknesses in order to make their surface as even as possible.
The proposed method was finally tested on a sample made of a canvas industrially prepared with a gypsum-based primer and covered with a layer of Cobalt Blue acrylic paint.  (Table 1), the five areas with the increasing thickness (1→5)are delimited by the horizontal dotted lines.

Micro-Raman Spectroscopy
In order to enhance the detection of Raman signals emitted by high fluorescent paints such as Cadmium Yellow (CY), Cadmium Red (CR), and Primary Red Magenta (PRM), Raman analysis was carried out with two devices using different excitation wavelengths.  (Table 1), the five areas with the increasing thickness (1→5)are delimited by the horizontal dotted lines.
The proposed method was finally tested on a sample made of a canvas industrially prepared with a gypsum-based primer and covered with a layer of Cobalt Blue acrylic paint.

Micro-Raman Spectroscopy
In order to enhance the detection of Raman signals emitted by high fluorescent paints such as Cadmium Yellow (CY), Cadmium Red (CR), and Primary Red Magenta (PRM), Raman analysis was carried out with two devices using different excitation wavelengths.
The first device is a Renishaw InVia 0310-02 System based on a continuous Nd:YAG laser excitation source at 532 nm. The diameter of the laser spot on the sample was diffraction limited to 1 µm by the objective lens (50×). The system is equipped with a Leica microscope (DMI 3000 M) and an electrically cooled CCD camera. The laser power was set between 0.15 and 0.30 mW.
The second device is a Renishaw Raman microscope RM1000 system coupled with an optical Leica DM LM microscope. The system is equipped with a refrigerated CCD camera and a CW He-Ne laser emitting at 632 nm, operating at a power of 3-30 mW, with a probing depth of 2 µm.
For all measurements, the spectral resolution was set at 4 cm −1 with an integration time in the range 10-60 s, the final spectra resulting from the accumulation of three individual ones.

Laser-Induced Breakdown Spectroscopy (LIBS)
The setup for LIBS analysis comprises a Q-switched Nd:YAG laser operating at its 4th harmonic at 266 nm, with a pulse duration of 15 ns and a repetition rate of 1 Hz. An f = 10 cm lens allows for fluences up to 6.6 J/cm 2 . The luminous emission was collected and dispersed by a 0.30 m spectrograph with a 1200 lines/mm grating (TMc300, Bentham, London, UK) coupled to an intensified charged coupled device (ICCD, 2151 Andor Technologies, Belfast, UK). Spectra were recorded at a 0.2 nm resolution, with a temporal gate of 3 µs and a delay of 500 ns, in order to avoid the continuous Bremsstrahlung emission. The laser beam was directed toward the surface of the samples at an angle of 45 • by using different mirrors. A cut-off filter of 300 nm was placed at the entrance of the spectrograph to reduce the scattered laser light and avoid second-order diffraction. The shot-to-shot fluctuation of laser pulse energy was less than 10%. The spectra resulted from summing the emissions of the ablation products after five successive laser pulses, a number that provided a good signal-to-noise ratio.

Spectral-Domain Optical Coherence Tomography (Sd-OCT)
The OCT device used in this study is a Thorlabs Telesto-II, working in the 1300 nm regime, with an axial resolution of 5.5 µm in air and a lateral resolution of 13 µm. The maximum field of view (FOV) is 10 × 10 mm 2 , with a 3.5 mm imaging depth. The system is controlled via a 64-bit software running on a high-performance computer. The 3D scanning probe with an integrated video camera allows for high-speed imaging (76 kHz) for rapid volume acquisition and live display. The sample stage provides XY translation and rotation of the sample along with the axial travel of the probe. The 2D tomograms were acquired with a pixel size of 3.5 µm 2 .

Fibre Optics Reflectance Spectroscopy (FORS)
FORS spectra were acquired with a Zeiss Multi-Channel Spectrometer, including two modules, the MCS601 UV-NIR C and the MCS611 NIR 2,2, with spectral sensitivity in the 190-1015 nm and 900-2200 nm range, and a spectral resolution of 0.8 nm/pixel and 5 nm/pixel, respectively. The size of the illumination spot was ∅ = 3 mm, and the illumination/observation geometry was 45 • /0 • . The output signal was processed through dedicated software, providing also CIEL*a*b* coordinates with standard D65 illuminant and 2 • observer. The calibration procedure was performed following CIE indications for non-contact spectrophotometric measurements by measuring a certified white 100% reflectance reference standard (Spectralon, Labsphere TM , North Sutton, PA, USA) and background noise.
For this specific application, it was essential to obtain spectra with reproducible reflectance, as the absolute intensity value is the core physical quantity to prove our thesis. Therefore, we defined an optimized procedure to both stabilize the signal intensity at each measurement and minimize fluctuations between the two modules. For each area (each of the five thicknesses), five reflectance spectra were measured in five different points with three consecutive acquisitions interleaved with dark correction by automatic closure of the device shutter.
For each area, five averaged spectra (over three acquisitions to check the signal stability) were obtained, which were then averaged to obtain a single spectrum for each thickness, together with the standard deviation. The reflectance of the white and black backgrounds was also measured following the same procedure. Therefore, a total of 1530 spectra were acquired on the painting mock-up (510 spectra for each measured point). Spectral data were processed in OriginLab environment.

Reflectance Imaging Spectroscopy (RIS)
The multispectral scanner [36] developed at CNR-INO operates in the range 395-2550 nm providing simultaneously 32 narrow-band images (16 VIS + 16 NIR) and pointwise spectral information. The simultaneous movement of both the lighting system and collecting optics, placed in a 45 • /0 • illumination/detection geometry, allows for uniform illumination and minimal heating of the surface. During the scanning, an autofocus system ensures the optimal target-lens distance. The system has a 250 µm sampling step (4 points/mm) and 500 mm/s maximum speed. A proper calibration procedure was performed by measuring a certified white 100% reflectance reference standard (Spectralon™) and background noise, following CIE indications for non-contact spectrophotometric measurements. The spectra reported in this work were extracted from the spectral cube using an in-house developed software. Each spectrum is averaged over an area of ∅ = 8 mm selected on each painted region. Given that the system has a 250 µm sampling step (4 points/mm), each spectrum is averaged over~805 pixels.

Conclusions
The chemical characterization of the considered ten acrylic paint materials was obtained by LIBS and micro-Raman spectroscopies, providing complementary information for the identification of elemental and molecular compounds in each paint and highlighting, in some cases, discrepancies with the chemical composition declared by the manufacturer.
The results obtained with the proposed method show that there is a clear correlation between the spectral response of the diffuse reflectance and the micrometric thickness of acrylic paint layers. Using two DRS modes (FORS and multispectral scanner), it was possible to identify significant spectral features (maxima and inflexion points) where the exponential decay of the reflectance on the layer thickness was evident. The resulting exponential fitting functions showed good agreement between the two techniques, both showing opposite trends depending on the substrate (white or black) for all paints. The fitting function defined for cobalt blue was used to estimate the thickness of a layer of paint spread on canvas, which was consistent with the thickness measured by OCT.
Measuring the thickness of a pictorial layer from its spectral characteristics can overcome some limitations of the OCT technique. Knowing a priori the refractive index of the material in order to correct optical distances in OCT tomograms is often not possible, especially for pigment mixtures typically found on paintings. Here, we have shown that the proposed method allows the measuring of the thickness of ten acrylic paints, including optically opaque pigments that cannot be measured by OCT, such as CY and CR (only in the presence of a black background), without the need to know their refractive index.
In this work, the feasibility of the method was demonstrated, assuming a single-layer homogeneous paint was present. For more complex situations, such as layers stratifications or pigment mixtures in painted artworks, more advanced techniques for feature extraction can be used [37]. Moreover, a systematic study based on specific spectral reflectance models simulating the painting composition and structure should be developed, considering the optical behaviour of each component. The presence of multilayers or pictorial mixtures, commonly found in paintings, requires the creation of a much larger database. In this view, approaches based on AI could be tested. An interesting development of the technique is the possibility of obtaining thickness maps on painted surfaces using hyperspectral reflectography, which offers spectral resolution comparable to FORS, but over large scanning areas.
Supplementary Materials: The following supporting information can be downloaded at: https: //www.mdpi.com/article/10.3390/molecules28124683/s1, Figure S1: Chemical characterization of the acrylic paints with LIBS and micro-Raman spectroscopy.

Data Availability Statement:
The data presented in this study are available on reasonable request from the corresponding author.